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We present Monte Carlo simulations of the spanning-forest model (q — * limit of the ferromagnetic 
Potts model) in spatial dimensions d — 3,4, 5. We show that, in contrast to the two-dimensional 
case, the model has a "ferromagnetic" second-order phase transition at a finite positive value w c . 
We present numerical estimates of w c and of the thermal and magnetic critical exponents. We 
conjecture that the upper critical dimension is 6. 
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The Potts model [l], Q plays an important role in the 
modern theory of phase transitions and critical phenom- 
ena, and is characterized by two parameters: the number 
q of Potts spin states, and the nearest-neighbor coupling 
v — e' 3,7 — 1. Initially q is a positive integer and v is a real 
number in the interval — 1 < v < +00, but the Fortuin- 
Kasteleyn (FK) representation 0] shows that the parti- 
tion function Zc{q,v) of the g-state Potts model on any 
finite graph G is in fact a polynomial in q and v. This 
allows us to interpret q and v as taking arbitrary real 
or even complex values, and to study the phase diagram 
of the Potts model in the real (q, w)-plane or in complex 
(q, u)-space. In particular, when q, v > the FK rep- 
resentation has positive weights and hence can be inter- 
preted probabilistically as a correlated bond-percolation 
model: the FK random- cluster model (J|. In this way we 
can study all positive values of q, integer or noninteger, 
within a unified framework. 

In two dimensions, the behavior of the ferromagnetic 
(v > 0) Potts/random-cluster model is fairly well un- 
derstood, thanks to a combination of exact solutions 5j, 
Coulomb-gas methods [6| and conformal field theory 7j. 
But in dimension d > 3, many important aspects remain 
unclear: the location of the crossover between second- 
order and first-order behavior [8| ; the nature of the criti- 
cal exponents and their dependence on q; the value of the 
upper critical dimension for noninteger q; and the quali- 
tative behavior of the critical curve v c {q) near q = 0. 

Interesting special cases of the random-cluster model 
arise in the limit q — > 0. In particular, the limit q, v —y 
with w = v/q held fixed gives rise to a model of spanning 
forests, i.e. spanning subgraphs without cycles, in which 
each occupied edge gets a weight w Very recently, 
it was shown — generalizing Kirchhoff 's matrix-tree 
theorem 11[ — that this spanning-forest model can be 



mapped onto a fermionic (Grassmann) theory involving a 
quadratic (Gaussian) term and a special nearest-neighbor 
four-fermion term. Moreover, this fermionic model pos- 
sesses an OSP(l\2) supersymmetry and can be mapped, 
to all orders of the perturbation theory in powers of 1/w, 



onto an TV- vector model (O(iV)-invariant a- model] ana- 
lytically continued to N = — 1. It follows that, in two 
dimensions, the spanning-forest model is perturbatively 
asymptotically free, in close analogy to (large classes 
of) two-dimensional er-models and four-dimensional non- 
abelian gauge theories. In particular, the only ferromag- 
netic (w > 0) critical point lies at w c = +00, in agree- 
ment [121 ] with the exact solutions on the square, triangu- 
lar and hexagonal lattices [5| showing that v c (q) oc q 1 ^ 2 
as q i 0. 

In this Letter we study the spanning-forest model in 
spatial dimensions d > 3, using Monte Carlo methods. 
We will show that, in contrast to the two-dimensional 
case, the model has a "ferromagnetic" second-order phase 
transition at a finite positive value w c , and we will esti- 
mate the thermal and magnetic critical exponents as well 
as a universal amplitude ratio. It follows that v c (q) oc q as 
q i 0. Indeed, we see the present study of the spanning- 
forest model as the first step in a comprehensive study of 
the random-cluster model as a function of (noninteger) q. 

For the random-cluster model with q > 1, a collective- 
mode Monte Carlo algorithm has recently been invented 
by Chayes and Machta [13] ; it generalizes the well-known 
Swendsen-Wang algorithm [14| and reduces to (a slight 
variant of) it when q is an integer. But for q < 1, the 
only available algorithm seems to be the Sweeny algo- 
rithm [15| . which is a local bond- update algorithm. Or- 
dinarily one would expect such a local algorithm to ex- 
hibit severe critical slowing-down, at least when the spe- 
cific heat is divergent [l6[ • But the random-cluster model 
with q < qo(d) s=s 2 has a non-divergent specific heat (i.e., 
critical exponent a < 0), which suggests that the critical 
slowing-down might not be so severe after all. Indeed, 
our numerical studies of the spanning-forest model (i.e., 
the q — y limit) in dimensions d = 2,3,4,5 strongly 
suggest that there is no critical slowing-down, i.e., the 
dynamic critical exponent z cxp associated to the expo- 
nential autocorrelation time is zero. Better yet, the ex- 
ponent £int,© associated to the integrated autocorrelation 
time [TtJ turns out to be negative for "global" observables 
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such as the mean-square cluster size; that is, one "effec- 
tively independent" sample can be obtained in a time 
much less than a single "sweep" — a kind of "critical 
speeding-up" . 

On the other hand, the Sweeny algorithm for q =/= 1 re- 
quires a non-local connectivity check each time one tries 
to update a single bond. If done in the naive way (e.g., by 
depth- first or breadth-first search), this would require a 
CPU time of order the mean cluster size \ ^ = L~ 2 
per "hit" of a single bond, leading to a severe "compu- 
tational critical slowing-down" . Recent work by com- 
puter scientists on dynamic connectivity algorithms fli| 
shows how this can be reduced to (logL) p , but at the 
expense of fairly complicated algorithms and data struc- 
tures. We therefore adopted an intermediate solution: 
a simple "homemade" dynamic connectivity algorithm 
that empirically has a slowing-down L~ - 7 . The details 
of this algorithm, along with measurements of the dy- 
namic critical behavior of the Sweeny algorithm in the 
spanning-forest limit, will be reported separately flit ]. 

We simulated the spanning-forest model in dimensions 
d = 3, 4, 5 on hypercubic lattices of size L d with periodic 
boundary conditions. We measured the cluster-size mo- 
ments S k = Eciustorsc #( c ) fc for & = 0, 2, 4. We focussed 
attention on the ratio R = (54)7(51), which tends in the 
infinite- volume limit to in a disordered phase and to 1 
in an ordered phase, and is therefore diagnostic of a phase 
transition. We also studied (52 ) in order to estimate the 
magnetic critical exponent. 

In each dimension, we began by making a "coarse" set 
of runs covering a wide range of w values, using modest- 
sized lattices and modest statistics. If the plots of R ver- 
sus w indicated a likely phase transition, we then made 
a "fine" set of runs covering a small neighborhood of the 
estimated critical point, using larger lattices and larger 
statistics. Finally, using the results from these latter 
runs, we made a "super-fine" set of runs extremely close 
to the estimated critical point, using as large lattices and 
statistics as we could manage, with the goal of obtaining 
precise quantitative estimates of the critical point w c and 
the critical exponents. The complete set of runs reported 
in this Letter used approximately 7 years CPU time on 
a 3.2 GHz Xeon EM64T processor. 

The "coarse" plot of R versus w for dimension d = 3 
and lattice sizes 6 < L < 32 is shown in Figure [TJ and 
shows a clear order-disorder transition at w c w 0.43. The 
corresponding "super-fine" plot, for lattice sizes 32 < 
L < 120, is shown in Figured We fit the data to Ansatze 
obtained from 

R = R c + cn(w - w c )L 1/u + a 2 (w - w c ) 2 L 2/v 

+ b x L- Ui + b 2 L-^ + ... (1) 

by omitting various subsets of terms, and we systemat- 
ically varied L m in (the smallest L value included in the 
fit). We also made analogous fits for (S^i/Ipl". Com- 
paring all these fits, we estimate the critical point w c = 
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Figure 1: Coarse plot of R versus w for spanning forests in 
dimension d = 3 and lattice sizes 6 < L < 32. 
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Figure 2: "Super-fine" plot of R versus w for spanning forests 
in dimension d — 3 and lattice sizes 32 < L < 120. 

0.43365 ± 0.00002, the critical exponents v = 1.28 ± 0.04 
and j/v — 2.1675 ± 0.0010, and the universal amplitude 
ratio R c = 0.8598 ± 0.0003 (68% subjective confidence 
intervals, including both statistical error and estimated 
systematic error due to unincluded corrections to scal- 
ing). A finite-size-scaling plot using these parameters is 
shown in Figure [31 A "coarse" plot of (5 2 )/L 7/l/ using 
the estimated value of j/v is shown in Figure |U 

The "coarse" plots of R versus w for dimensions d = 
4, 5 are shown in Figures [5] and [U respectively. Once 
again they show a clear order-disorder transition. For 
lack of space, we refrain from showing the corresponding 
"super-fine" plots (which use lattice sizes up to 64 4 and 
20 5 ) and simply give the results of fits to Ansatze of the 
general type ([1]). In dimension d — 4, we estimate w c = 
0.210302 ± 0.000010, v = 0.80 ± 0.01, -y/v = 2.1603 ± 
0.0010 and R c = 0.73907 ± 0.00010. In dimension d = 5, 
we estimate w c = 0.14036 ± 0.00002, v = 0.59 ± 0.02, 
7/f = 2.08 ± 0.02 and R c = 0.625 ± 0.015. 
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Figure 3: Finite-size-scaling plot of R versus (w — w c )L 1 ^ I/ , Figure 6: Coarse plot of R versus w for spanning forests in 
with w c = 0.43365 and v = 1.28, for spanning forests in dimension d = 5 and lattice sizes 4 < L < 12. 
dimension d — 3 and lattice sizes 8 < L < 120. 




0.2 0.3 0.4 0.5 



w 

Figure 4: Plot of (S 2 )/L y/,/ versus w, with j/v = 2.1675, for 
spanning forests in dimension d = 3 and lattice sizes 6 < L < 
32. 
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Figure 5: Coarse plot of R versus w for spanning forests in 
dimension d — 4 and lattice sizes 4 < L < 20. 
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q = 2 


d = 2 


f = oo 
j/u = 2 


^ = 4/3 
7/1/ = 43/24 


v = \ 

f/v = 7/4 


d= 3 


k = 1.28(4) 
7/i/ = 2.1675(10) 


i/ = 0.874(2) 
-y/v = 2.0455(6) 


v = 0.6301(5) 
7/1/ = 1.9634(5) 


d= 4 


v = 0.80(1) 
j/v = 2.1603(10) 


v = 0.689(10) 
7/i/ = 2.094(3) 


f = 1/2 (log) 
7/v = 2 (log) 


d = 5 


v = 0.59(2) 
7/v = 2.08(2) 


v = 0.57(1) 
7/^ = 2.08(2) 


v = 1/2 
7/^ = 2 



Table I: Critical exponents v and 7/1/ versus q and d. d = 2: 
presumed exact values [2(j. d = 3,4,5, g = 0: this work, 
d = 3, q = 1: J2I d = 3, 9 = 2: (22(|. d = 4, g = 1: [H, H3]. 
d = 5, g = 1: 24, 25]. d — 4,5, q — 2: presumed exact values. 



In Table U we summarize the estimated critical ex- 
ponents for ferromagnetic Potts models with q = 
(this work), 1 (percolation) and 2 (Ising) in dimensions 
d = 2, 3, 4, 5. It is evident that v varies quite sharply as 
a function of q and d, while 7/1/ varies much more slowly. 
The dimension-dependences of v and 7/1/ for q = are 
consistent with the conjecture that they are tending to 
the mean-field values 1/2 and 2 in dimension d = 6, just 
as they do for q = 1. This in turn supports the more 
general conjecture that the upper critical dimension is 6 
for all random-cluster models with < q < 2, and is 4 
only when q = 2. 

This conjecture is supported by a field-theoretic 
renormalization-group calculation in dimension d = 6 — e 
through order e 3 [261 ] in which q — 2 plays a distinguished 
role (all the correction terms vanish there). Specializing 
to q = 0, we have 

7/^2 + ± + ^-f^-^V + ^) 
" 15 225 V 625 16875/ V ' 
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2 + 0.066667e + 0.031111e 2 - 0.034065e 3 + 0(e 4 ) 



1/u 



2-1-1- 
3 30 



173 \ 



4C(3) 



125 27000 



+ 0(e 4 ) 



= 2 - 0.333333e - 0.033333e 2 + 0.032058e 3 + 0(e 4 ) 

These series seem rather difficult to resum, especially 
when e > 2, but they are in qualitative agreement with 
the exponents listed in Table fl] Moreover, a slightly bet- 
ter agreement can be obtained by imposing the known 
exact values at e — 4 on the interpolating function. 

Details of these simulations and their data analysis, in- 
cluding analysis of universal amplitude ratios other than 
R, will be reported separately [lj|. 
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